Magnetism in a Ti02/LaA103 heterostructure: an ab initio study about the role of 

oxygen vacancies 
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In this work we study the electronic structure and magnetism of a Ti02 film grown on another 
non-magnetic oxide such as a LaAlOa (001) substrate, concentrating on the role played by structural 
relaxation and oxygen vacancies. Using Density Functional Theory ab-initio methods, we study the 
free-standing anatase film as well as the interfaces with either the LaO or AIO2 planes of LaAlOa, 
focusing on the possibility of magnetic solutions. Our results show that the interface LaO/Ti02 
is favored against the A102/Ti02 one if no oxygen vacancies are present in the interface whereas 
the contrary happens when there are oxygen vacancies. In both cases, the cohesive energy is of the 
same order of magnitude but only the A102/Ti02 interface presents an stable magnetic solution. 

PACS numbers: 73.20.-r, 75.70.-i 
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I. INTRODUCTION 

The development of modern microelectronics de- 
vices requires the knowledge of the electronic struc- 
ture near interfaces. In this respect, oxide materials 
combine many of the important properties of semicon- 
ductors and add novel phases such as metal-insulator 
transitions^, superconductivitj*^, quantum Hall effecfci, 
orbital reconstruction-, colossal magnetoresistance^,etc. 
One interesting phenomena reported recently is the 
change in the properties of complex oxides when they 
are arranged in a heterostructure. For example, LaAlOa 
(LAO) and SrTiOa (STO) are both insulators and non- 
magnetic in their bulk form but one possible interface 
between them is conductive^ and ferromagnetio^. The 
origin of the ferromagnetism and the conductivity is not 
known, it has been attributed to the polar discontinu- 
ity between the charged layer (LaO) in contact with the 
neutral layer (Ti02), that would produce a transfer of 
electrons at the interface, and also to the presence of 
oxygen vacancies created during film growt h^'^d^ . 

The previous findings may be related to the recent 
intense search for ferromagnetism above room temper- 
ature in dilute magnetic semiconductors and insulators. 
The fabrication of these materials offers exciting possi- 
bilities for spintronic devices and the use of wide gap ox- 
ides is appealing for magnetooptical devices. In that di- 
rection, room-temperature ferromagnetism has been ob- 
served in both dope d^^i^^ and undoped insulating oxide 
thin films such as Ti02, ZnO and HfO a^'^i^^'^^ grown over 
other oxides or semiconductors. This magnetic order can 
be weak and the determination of its intrinsic charac- 
ter may require complementary, independent and care- 
ful techniques to be probe d^-'^d^d'^ . Many authors have 
attributed an important role to oxygen vacanciesiS. and 
other defectsi^i22, to produce this magnetic ordering, but 
the effect of the substrate or the interface has not been 
taken into account in the interpretations. 

Considering the open questions in the aforementioned 



systems, in this work we explore the possibility of observ- 
ing magnetism in an heterostructure formed by a simple 
oxide such a Ti02 anatase and a LAO substrate. Exper- 
iments show that large terraces are formed on the sub- 
strate surface, and that the growth is epitaxial in the 
(001) direction2i. However, it is nor clear if the surface 
is mostly composed of LaO planes or of AIO2 planes^i^. 
XPS experiments show that there are Ti+^ and Ti+^ ions 
in these films, which may be due to the presence of oxygen 
vacancies^ or to the polar catastrophe. The magnitude 
of the magnetization has been found to be proportional 
the vacancy concentration in some experiments, for ex- 
ample in ref." We therefore concentrate on the effect of 
oxygen vacancies in the electronic structure and on the 
interfacial structural relaxation, as it is expected that 
the electric fields due to the dipole layer at the inter- 
face may lead to important changes in the positions of 
the atoms with respect to the bulk materials, as it hap- 
pens in other interfaces^. For this purpose we perform 
ab-initio calculations for the different interfaces present 
in the Ti02/LA0 heterostructure using computational 
codes such as SIESTA^^ or Wien2k2& and also different 
unit cells. 

We first study separately the component systems: 
Ti02 and LAO in bulk and as free standing slabs, rep- 
resenting the deposited anatase film and the substrate, 
respectively. As the ab-initio codes allow only the study 
of three dimensional periodic systems, the films are repre- 
sented by repeated slabs with enough free space between 
them. For the heterostructure, two different types of unit 
cell were used: trilayer slabs LA0/Ti02/LA0 and super- 
lattices ...Ti02/LAO/Ti02/LAO... 

The paper is organized as follows: In Sec. [TT]we very 
briefly outline our calculational approach. In Sec lIIII we 
discuss the bulk and surface electronic properties of the 
component materials. In Sec. IIVI we present the results 
for the trilayers and superlattices and in Sec. IVlour con- 
clusions. 
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fin tin radii used were R^i^ 



bohr, RfJ 



2.5 bohr, R*7 



II. CALCULATION DETAILS 

The calculations are performed within the density 
functional theory (DFT)21 using the full potential aug- 
mented plane waves method, as implemented in the 
Wien2k code^^. We use the local density approxima- 
tion (LDA)22. for the exchange and correlation, and small 
muffin tin radii to allow for lattice relaxation. The muf- 

= 1.7 bohr, R^^ = 1.4 
: 1.7 bohr. The number of 
plane waves in the interstitial region is mostly restricted 
to RKMax=6 although in many cases it was increased to 
RKMax=7 and no qualitative difference was observed. 
The effect of increasing the number of k-points was also 
checked. In some cases we relax the structure with the 
Siesta code^^ as it allows for further cell deformation. 



III. TIO2 ANATASE AND LAO: BULK AND 
SURFACE PROPERTIES 

Ti02 anatase and LAO are both conventional band in- 
sulators, with in-plane bulk lattice parameter a,—3.79A^ 
(Fig. [1]). The small mismatch between these in-plane 
lattice constants allows for a good epitaxial experimen- 
tal outcome^. 



gap and the vacancy levels very close to the conduction 
band are typical features of the LDA approximation. 
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FIG. 1: Bulk unit cells of the two component systems. The 
arrows on the oxygen atoms in the anatase structure indicate 
that they are ofl-plane with respect to the Ti atoms. 
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FIG. 2: Density of States of: (a) bulk Ti02 anatase (b) 
bulk LAO, (c) and (d) Superlattices either with AIO2 or LaO 
planes in the interfaces (relaxed structures). Fermi energy is 
set to O (in all DOS graphs). 
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Ti02 anatase presents a 3.2 cV band gap between filled 
oxygen 2p bands and unfilled Ti 3d conduction bands . 
The calculated band gap is 2.2 eV as shown in Fig. [H^a). 
When oxygen vacancies are formed in the bulk anatase 
structure the calculations show that the system becomes 
metallic and a vacancy level appears inside the gap, very 
close to the conduction band (Fig. [3]). The small band 



FIG. 3: Density of states of bulk anatase with one oxygen 
vacancy (Tiie O31). Notice the vacancy energy level very 
close to the conduction band. 

LAO presents a 5.6 eV gap between filled oxygen 2p 
bands hybridized with Al p states and unfilled conduction 
bands composed mostly of La 4/ states. The calculated 
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gap is 4 eV, as shown in Fig. [2tb), and when oxygen va- 
cancies are created in this system the calculated vacancy 
levels appear well inside the band gap. 

Summing up, both materials are non-magnetic insu- 
lating oxides in bulk and if oxygen vacancies are present 
they keep their non-magnetic character although the elec- 
tronic structure changes, as vacancy levels appear inside 
the band gap. 

In order to get some insight about how these bulk prop- 
erties change in the heterostructure, we study slabs of 
the component materials along the (001) direction. The 
anatase slab with 5 or 9 Ti02 layers is used to model 
the properties of the free surface of an anatase film. A 
different DFT calculation, using the SIESTA code, al- 
lows the unit cell to change its shape under relaxation, 
which is not the case with the Wien2k code. This results 
in a strong deformation of the (001) surface, that goes 
from a square into a rectangular shape, and also makes 
the slab thinner. The free standing anatase film would 
therefore deform in that way, but as we are studying its 
epitaxial growth over LAO we keep the in plane lattice 
constant fixed and the cell shape square. Relaxation is 
only considered in the z direction, perpendicular to the 
slab. With this constraint we create oxygen vacancies in 
the surface of the slab and relax the atomic positions us- 
ing the Wien2k code. The vacancies seem to be crucial 
for the appearance of magnetism: on one hand, mag- 
netic solutions appear if there are oxygen vacancies at 
the surface and on the other hand the value of the mag- 
netic moment increases with the number of vacancies. In 
our calculations, one vacancy per surface (50%) gives af- 
ter relaxation a magnetic moment of 0.3 /zb, localized 
at the Ti atom in that surface, while two vacancies per 
surface give 2 fiB- In Fig- El we show the charge and spin 
density maps near the surface for the anatase slab with 
two vacancies per surface (all the surface oxygen atoms 
removed). The system is clearly magnetic, and the mag- 
netic moment is mostly localized in the Ti superficial 
atom. The effect of structural relaxation is to move this 
Ti atom closer to the remaining oxygen neighbor, which 
is apical and in the subsurface layer, thus making the 
slab thinner. In fact, this distance decreases from 1.97 
A in bulk Ti02 to 1.89 A for a surface with no oxygen 
vacancies and to 1.75 A if there are oxygen vacancies in 
the surface layer as in Fig. ID In the slab with 5 Ti02 
layers the two surfaces interact in such a way that the 
magnetic moments are smaller than in the 9-layer slab. 

The other building block of the heterostructures is the 
slab of LAO. Along the (001) direction this slab alter- 
nates layers of AIO2 (negatively charged) and LaO (pos- 
itively charged). Therefore, an integer number of unit 
cells in the z direction would produce a large dipole and 
it is more convenient to use an odd number of layers. 
For the superlattices we consider 5 layer slabs with either 
LaO or AIO2 termination. In both cases the system be- 
comes metallic within the LDA approximation and the 
interesting change with respect to the bulk density of 
states, shown in Fig. [H^b), is that the LaO termination 
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FIG. 4: Charge and spin density of the surface of a Ti02 slab 
with no oxygen atoms at the surface. Notice that only the 
surface Ti atom has a considerable magnetic moment. 



shifts the unoccupied surface La states substantially to- 
wards the Fermi energy. This can give as a result a large 
hybridization with the unoccupied Ti states of Ti02, but 
will not modify the occupied states. 



IV. TIO2/LAO SUPERLATTICES 

The system we propose to study is a Ti02 (001) film 
with anatase structure grown epitaxially over a LAO 
(001) substrate. For this purpose we perform calcula- 
tions on the superlattices shown schematically in Fig. [5] 
We have considered the 2 cases: the Ti02 plane facing 
the AIO2 or the LaO termination. In both situations we 
assume that the oxygen atoms of the LAO surface face 
the Ti atoms of the anatase surface. The unit cell we use 
consists of 5 layers of each material and contains 27 or 
28 atoms. The lateral lattice parameter is kept fixed at 
the experimental value for LAO but the distance between 
the two slabs composing the heterostructure is obtained 
by minimizing the total energy. The fact that the total 
energy has a minimum when changing this distance veri- 
fies that there is bonding between the two materials and 
also gives the optimal value for that distance. This fixes 
the size of the superlattice unit cell along the z axis for 
each case, for example we notice that it changes when 
there are oxygen vacancies. Afterwards, all the atoms 
are allowed to relax inside the unit cell until the force on 
each one is less than 2 me v/A. 

We find that there is a competition between the ten- 
dency of the anatase film to shrink, as in the free standing 
slab, and the tendency to bond to LAO. In particular, if 
an oxygen from the anatase surface faces an Al atom the 
bond is strong. Also, there is an electrostatic effect due 
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FIG. 5: Schematic diagrams for the unit cells of the superlat- 
tices studied. Notice that there are two types of A102/Ti02 
interfaces, due to the off-plane oxygen atoms in the anatase 
structure (Fig. [T]), but only one type of LaO/Ti02 interface. 



to the charge in each layer of LAO, which is of course 
different for a neutral Ti02 layer and for a charged layer, 
when oxygen vacancies are created in the anatase surface. 

Fig- IHc) and (d) show the density of states of the 
two superlattices without vacancies and after relaxation. 
These graphs give some insight about electron transfer 
or the possibility of oxygen vacancy formation, in agree- 
ment with experiments^. In the case of the A102/Ti02 
interface, the states readily available near Ep are 0-2p 
states. This interface will then be prone to the formation 
of oxygen vacancies in order to allow for a compensation 
of the polar discontinuity present at the interface. For 
the LaO/Ti02 interface instead, the states near E p are 
mainly comprised of Ti 3ci-states. This accessibility to a 
mixed valency of Ti allows for a transference of electrons 
across the interface to overcome the electrostatic poten- 
tial divergence. Both interfaces are metallic without va- 
cancies. This metallicity is not altered by the formation 
of oxygen vacancies or by structural relaxation. 



A. A102/Ti02 heterointerface 
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FIG. 6: Charge density for the superlattice with A102/Ti02 
interfaces and no vacancies. The two interfaces are differ- 
ent: one can see that the distance from the surface Al to the 
neighbor anatase oxygen is shorter in interface 1. 



ficial atoms and do not disappear, they only decrease 
somewhat. One particular example of the structural re- 
laxation is shown in Fig. [71 indicating schematically how 
each of the atoms moves in the magnetic interface. The 
atom that has a larger displacement is the anatase oxygen 
facing Al, it bonds more strongly with Al than with Ti. 
Relaxation changes the buckling of the anatase surface. 



Due to the anatase structure there are two types of 
A102/Ti02 interfaces. This can be inferred from Fig. [1] 
and Fig. [D^a), but it is shown more clearly in Fig. [6l To 
simulate growth conditions in a typical experiment*^, we 
introduced oxygen vacancies at both layers of the inter- 
face and found that it is more favorable to perform them 
in the Ti02 layer rather than in the AIO2 plane. 

When the heterostructure is formed, the interface be- 
comes metallic; with no oxygen vacancies the Fermi level 
lies in the valence band (Fig. [Ifc)) and with vacancies it 
moves to the conduction band. 

We find magnetic solutions in the unrelaxed geome- 
try, with and without vacancies, but the magnetic mo- 
ments decrease with structural relaxation as can be seen 
in Table 1. Without vacancies the moments are located 
in some of the oxygen atoms and disappear with relax- 
ation. With vacancies they are located in the Ti super- 
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FIG. 7: Schematic diagram showing the atomic relaxation on 
the magnetic interface (namely, interface 2 in Fig. O when 
there is one oxygen vacancy at the anatase surface. The move- 
ment of each atom (in A) is indicated when it is significant 
(more than 0.05 A) 

Fig. [5] shows the charge and spin densities of the su- 
perlattice, when there are no oxygen atoms in the Ti02 
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side of the interfaces (2 vacancies). Although both inter- 
faces have the same number of missing oxygen atoms and 
first neighbours, only one of them is magnetic, depending 
on the relative orientation of the second neighbor layers. 
This points to the new result that magnetism is not only 
due to interactions between interfacial atoms but also 
further neighbours are important. 

The corresponding density of states is shown in Fig. [9] 
where the Ti 3c? states are marked and they are mainly 
composed of dx2-y2 and d^j, states. 
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FIG. 9: Density of states for the case of Fig. |8] with the 
relaxed structure. The states corresponding to the Ti at in- 
terface 2 are highlighted. 



lency state for the Ti atoms at the interface without va- 
cancies while in the previous superlattice it was due to 
the vacancies. 

Atomic relaxation in this interface is considerably 
smaller than in the previous one. 



FIG. 8: Charge and spin density maps for the superlattice 
with AI02/Ti02 interfaces and 2 vacancies. The effect of 
relaxation is too small to be seen in this scale. Notice that 
only one of the interfaces is magnetic. 



B. LaO/Ti02 heterointerface 

To study this heterointerface we use the supercell de- 
picted in Fig. El^b) and perform oxygen vacancies to 
simulate experimental growth conditions, ft is more fa- 
vorable to perform them in the Ti02 layer compared to 
the LaO face. Introducing vacancies in this case is over- 
all less favourable compared to the other interface and 
magnetic solutions are only found with many oxygen va- 
cancies present, as shown in Table 1. As in this case 
both interfaces of the superlattice are identical, the mag- 
netic moment per Ti atom at the surface is considerably 
smaller than in the previous interface. In Fig. ^d) we 
can see that hybridization of La 4/ with Ti 3d electrons is 
large, probably exaggerated by the LDA approximation, 
as was mentioned in Refi^. Comparing the valence elec- 
tron charge inside the muffin-tin spheres for Ti with those 
of their counterparts in bulk, we corroborate a mixed va- 



C. Cohesive energies 

As a very first approximation one could estimate the 
cohesive energies of the systems in their relaxed struc- 
tures by subtracting the energy of the conforming atoms 
from the total energy of the system. By doing so wc ob- 
tain that a 5 layer slab of LAO prefers the termination 
AIO2 over LaO. The cohesive energy of a superlattice 
with A102/Ti02 interfaces is also larger than that with 
LaO/Ti02 interfaces. 

A better approximation to the cohesive energy of a 
superlattice, given in Table I, is to subtract from its to- 
tal energy the energy of the relaxed slabs of Ti02 and 
LAO that compose it. This gives a more interesting re- 
sult: the supercell with LaO/Ti02 interface has more 
cohesive energy if no vacancies are present while the one 
with A102/Ti02 interface is preferred if there are oxygen 
vacancies. This last case presents a magnetic moment, lo- 
calized at the superficial Ti atom. It is important to note 
that the energy difference between cases with one and 
two vacancies is order of magnitude smaller compared to 
the cases with no vacancies. Therefore, with our results 
we can draw qualitative conclusions comparing cases ei- 
ther with or without vacancies but not referring to the 
concentration of vacancies. 
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case 


11/ Al 


T'; /T o 




E (cV) M (ms) 


E(cV) /X (ms) 


no vac 


4.89 0.75 ^ 0.01 


6.86 


1 vac 


6.65 0.38 ^ 0.26 


6.01 


2 vac 


6.62 1.36^ 1.28 


5.94 1.05-> 0.57 



TABLE I: Cohesive energies and total magnetic moments per 
unit cell for the two types of superlattices considered, with 
or without vacancies (vac). The effect of lattice relaxation 
is indicated by the arrows for the magnetic moments. The 
number of vacancies is the number of missing oxygen atoms 
from the anatase side of each interface. 



CONCLUSIONS 



terface predominates. Using energetic arguments we find 
that Ti02 deposition over botli LAO layers could occur 
with equal probability but with different oxygen content 
in each case. The interface with AIO2 has the largest 
cohesive energy when there are oxygen vacancies, while 
when there are no vacancies, the interface with LaO is 
preferred. Therefore, if the more abundant terraces have 
LaO termination, the magnetism observed experimen- 
tally would not be interfacial. It could be due to va- 
cancies in the free surface of the anatase film, and the 
role of the substrate would be only to fix the cell size. 
If on the contrary the more abundant interface would be 
with Al, and the system is grown under low oxygen pres- 
sure, there could be a magnetic moment localized at the 
interface. 



Our calculation shows that magnetism appears at the 
Ti02 anatase surface and at the Ti02/A102 interface 
when there are enough oxygen vacancies. Structural re- 
laxation is found to be significant and produces an im- 
portant effect on the superficial Ti magnetic moment, di- 
minishing it. The experimental substrate LAO, presents 
large terraces in its surface but it is not clear which in- 
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